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RESUMEN 

Favor de proporcionar un resumen en espanol. If you are unable to translate your abstract 
into Spanish, the editors will do it for you. Recent observation of 850/mi sub-mm polarization from T 
Tauri disks opens up the possibility of studying magnetic field structure within protostellar disks. The degree 
of polarization is around 3 % and the direction of polarization is perpendicular to the disk. Since thermal 
emission from dust grains dominates the spectral energy distribution at the sub-mm/FIR regime, dust grains 
arc thought to be the cause of the polarization. We discuss grain alignment by radiation and we explore the 
efficiency of dust alignment in T Tauri disks. Calculations show that dust grains located far away from the 
Central proto-star are more efficiently aligned. In the presence of a regular magnetic field, the aligned grains 
produce polarized emission in sub-mm/FIR wavelengths. The direction of polarization is perpendicular to the 
local magnetic field direction. When we use a recent T Tauri disk model and take a Mathis-Rumpl-Nordsieck- 
type distribution with maximum grain size of 500-1000 /um, the degree of polarization is around 2-3 % level 
at wavelengths larger than ~100/um. Our study indicates that multifrcquency infrared polarimetric studies 
of protostellar disks can provide good insights into the details of their magnetic structure. We also provide 
predictions for polarized emission for disks viewed at different wavelengths and viewing angles. 

ABSTRACT 

Recent observation of 850^tm sub-mm polarization from T Tauri disks opens up the possibility of studying 
magnetic field structure within protostellar disks. The degree of polarization is around 3 % and the direction 
of polarization is perpendicular to the disk. Since thermal emission from dust grains dominates the spectral 
energy distribution at the sub-mm/FIR regime, dust grains are thought to be the cause of the polarization. 
We discuss grain alignment by radiation and we explore the efficiency of dust alignment in T Tauri disks. 
Calculations show that dust grains located far away from the Central proto-star are more efficiently aligned. 
In the presence of a regular magnetic field, the aligned grains produce polarized emission in sub-mm/FIR 
wavelengths. The direction of polarization is perpendicular to the local magnetic field direction. When we use 
a recent T Tauri disk model and take a Mathis-Rumpl-Nordsieck-type distribution with maximum grain size of 
500-1000 fim, the degree of polarization is around 2-3 % level at wavelengths larger than ~100^m. Our study 
indicates that multifrcquency infrared polarimetric studies of protostellar disks can provide good insights into 
the details of their magnetic structure. We also provide predictions for polarized emission for disks viewed at 
different wavelengths and viewing angles. 

Key Words: polarization — accretion disk — dust, extinction — Stars: Pre-main sequence 



1. INTRODUCTION 

Recently, Tamura ct al. (1999) first detected po- 
larized emission from T Tauri stars, low mass proto- 
stars. They interpreted the polarization (at ~ 3 % 
level) in terms of thermal emission from aligned dust 
grains. Magnetic field is an essential component for 
grain alignment. If grains are aligned with their long 
axes perpendicular to magnetic field, the resulting 
grain emission has polarization directed perpendicu- 
lar to the magnetic field. 



x Dcpt. of Astronomy and Space Science, Chungnam Na- 
tional Univ., Korea (jcho@cnu.ac.kr). 

2 Astronomy Dept., Univ. of Wisconsin, Madison, 
WI53706, USA (lazarian@astro.wisc.edu). 



The notion that the grains can be aligned in re- 
spect to magnetic field can be traced back to the dis- 
covery of star-light polarization by Hall (1949) and 
Hiltncr (1949), that arises from interstellar grains. 
Historically the theory of the grain alignment was de- 
veloping mostly to explain the interstellar polariza- 
tion, but grain alignment is a much wider spread phe- 
nomenon (see Lazarian 2007 for a review). Among 
the alignment mechanisms the one related to radia- 
tive torques (RTs) looks the most promising. Wc 
invoke it for our calculations below. 

The RTs make use of interaction of radiation with 
a grain to spin the grain up. The RT alignment was 
first discussed by Dolginov (1972) and Dolginov & 
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Mytrophanov (1976). However, quantitative studies 
were done only in 1990's. In their pioneering work, 
Draine & Weingartner (1996, 1997) demonstrated 
the efficiency of the RT alignment for a few arbitrary 
chosen irregular grains using numerical simulations. 
This work identified RTs as potentially the major 
agent for interstellar grain alignment. A successful 
analytical model of RTs was suggested by Lazarian & 
Hoang (2007). Cho & Lazarian (2005) demonstrated 
the rapid increase of radiative torque efficiency and 
showed that radiative alignment can naturally ex- 
plain decrease of the degree of polarization near the 
centers of pre-stellar cores. Large grains are known 
to be present in protostellar disk environments and 
this makes the RT alignment promising. 

Roughly speaking, the efficiency of grain align- 
ment by RTs depends on two factors - the intensity of 
radiation and the gaseous drag. The latter depends 
on gas pressure. Therefore, the ideal condition for 
grain alignment by RTs is strong radiation and low 
gas pressure. 

In order to calculate efficiency of grain alignment, 
we need to know radiation intensity, gas density, and 
temperature in T Tauri disks. Recently proposed hy- 
drostatic, radiative equilibrium passive disk model 
(Chiang & Goldrcich 1997; Chiang et al. 2001, here- 
after C01) fits observed SED from T Tauri stars very 
well and seems to be one of the most promising mod- 
els. Here, passive disk means that active accretion 
effect, which might be very important in the imme- 
diate vicinity of the central star, is not included in 
the model. In this paper we adopt the model in C01. 

In this paper, we briefly discuss polarized FIR 
emission arising from aligned dust grains by radiative 
torque in T Tauri disks. Detailed calculations and 
discussions can be found in Cho & Lazarian (2007). 
In §11, wc discuss grain alignment in T Tauri disks. 
In §111, we give theoretical estimates for degree of 
polarization. In §IV, we discuss observational impli- 
cations. Wc give summary in §V. 

2. THE DISK MODEL USED FOR THIS STUDY 
We assume that magnetic field is regular and 
toroidal (i.e. azimuthal). We use a T Tauri disk 
model in C01. Figure Q] schematically shows the 
model. The disk is in hydrostatic and radiative equi- 
librium and shows flaring. They considered a two- 
layered disk model. Dust grains in the surface layer 
are heated directly by the radiation from the central 
star and emit their heat more or less isotropically. 
Half of the dust thermal emission immediately es- 
capes and the other half enters into disk interior and 
heats dusts and gas there. They assume that the 
disk interior is isothermal. 



In our calculations, we use a grain model similar 
to that in C01. We use an MRN-type power-law dis- 
tribution of grain radii a between a m i n (=0.01 \im 
for both disk interior and surface layer) and a max 
(=1000 [im for disk interior and = lfim for disk sur- 
face layer) with a power index of -3.5: dN cx a~ 3 5 da. 
As in C01 we assume that grain composition varies 
with distance from the central star in both disk in- 
terior and surface layer. We assume that grains 
in the surface layer are made of silicate only when 
the distance r is less than 6 AU, and silicate cov- 
ered with water ice when r > 6AU. We do not 
use iron grains for the immediate vicinity of the 
star. We assume that grains in the disk interior 
are made of silicate when r < 0.8AU and ice-silicate 
for r > 0.8AU. The fractional thickness of the wa- 
ter ice mantle, Aa/a, is set to 0.4 for both disk 
surface and disk interior. Unlike C01, wc use the 
refractive index of astronomical silicate (Draine & 
Lee 1984; Draine 1985; Loar & Draine 1993; see 
also Weingartner & Draine 2001). We take optical 
constants of pure water ice from a NASA web site 
(ftp:// climatel .gsfc.nasa.gov/wiscombe I . 




Fig. 1. A schematic view of the disk model (see C01). 
The surface layer is hotter and heated by the star light. 
The disk interior is heated by re-processed light from the 
surface layers. We assume that the disk height, H, is 
4 time the disk scale height, h. From Cho & Lazarian 
(2007). 

3/2 

The column density of the disk is Yi$r Al j with 
So = 1000g/cm 2 . Here tau is distance measured in 
AU. The disk is geometrically flared and the height 
of the disk surface is set to 4 times the disk scale 
height h. The disk inner radius is 2i?* and the 
outer radius is 100AU. The central star has radius 
of i?* = 2.5Rsun and temperature of T* = 4000K. 
Temperature profile, flaring of disk, and other details 
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of the disk model are described in C01. 

3. GRAIN ALIGNMENT BY RADIATION 

3.1. Polarized FIR emission from aligned grains 

Theories predict that grain alignment happens in 
such a manner that grain's long axis is perpendicu- 
lar to local magnetic field direction (see Figure [2|) . 
Grains are usually cold and emit infrared (IR) radi- 
ation. When an elongated grain emit IR radiation, 
the one with electric vector parallel to the grain's 
long axis is stronger (left panel of Figure [2]) . There- 
fore the direction of polarization is parallel to the 
grain's long axis, or perpendicular to the magnetic 
field. 

Then what is the condition of grain alignment? 
Roughly speaking, fast rotation is a necessary condi- 
tion for grain alignment. Then what is the condition 
for fast rotation? Of course, when the intensity of 
radiation is strong, grains can rotate faster. Gases 
density is also an important factor. Since gaseous 
drag slows down grain's rotation, low gas density 
(more precisely, gas pressure) is a favorable condi- 
tion for fast rotation. (Note however that this is just 
a possible parameterzation for RTs. See Lazarian & 
Hoang 2007 and the last paragraph of this section 
for details.) 

3.2. Radiative torque for large grains 

For most of the ISM problems, dust grains are 
usually smaller than the wavelengths of interest. 
However, this is no longer true in T Tauri disks 
because we are dealing with grains as large as ~ 
1000/im. To understand grain alignment in T Tauri 
disks we need to know radiative torque for large 
grains. Motivated by Figure [H we assume that the 
radiative torque 
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Qr = 



0(1) 
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if A ~ a 
if A > a, 
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where a is the grain size, Urad the energy density 
of the incident radiation, and A the wavelength of 
the incident radiation. Note that Qr ~ 0(1) when 
A ~ a. 

3.3. Rotation rate of dust grains by radiative torque 

We assume that to align grains RTs should spin 
grains suprathermally. Detailed theory of grain 
alignment can be more complicated (see recent study 
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Fig. 3. Behavior of Torque. Torque is ~ O(l) when 
A ~ a, where a is the grain size. Roughly speaking, 
torque oc (A/a)~ 3 . The results for large grains. Data 
from Lazarian & Hoang (2007). 

by Lazarian & Hoang 2007, 2008; Hoang & Lazarian 
2007). 

After some modifications, equation (67) in 
Draine & Weingartner (1996) reads 

2 1n 9 a l ( u rad \ 2 ( A 

= 4 - 72xl0 V 3a - 5 U^J '- 
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Tdrag 
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where Qr = Qr ■ &i and ai is the principal axis 
with largest moment of inertia, hh is the hydrogen 
number density, u ra d is the energy density of the ra- 
diation field, S rs 2, ai « 1.745, p% = rho/3gcm~ 3 , 
0,5 = a/10~ 5 cm, and u>t is the thermal angular fre- 
quency, which is the rate at which the rotational 
kinetic energy of a grain is equal to kT/2. The 
timescales Td ra g,gas and T drag ^ m are the damping 
time for gas drag and for electromagnetic emission, 
respectively, and they satisfy the relation T^ rag — 

T drag,e m + T drag,gas ( SCC DrainC & Weingartner 1996 

for details). As we discussed in the previous sub- 
section, Qr is of order of unity when A ~ a and 
declines as (A/a) increases. From this observation, 
we can write 



l^rad 
UJT 



for A > a, where 



Urad 

U>rn. 



Or, 
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LOT 



4.72 x 10 9 ^p 3 a-5 , 
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Elongated grains emit polarized IR 




radiation 



Stronger 



condition for alignment = fast rotation 



Fig. 2. Grain alignment and polarized emission. Left Emission from an elongated grain. When electric field is parallel 
to the grain's long axis, the radiation is stronger. Right: Conditions for grain alignment. Roughly speaking, low gas 
density and strong rotation are favorable conditions for grain alignment. 
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The limitation of this approach is that the ampli- 
tude values of the RTs are used to parameterize the 
alignment. In fact, Lazarian & Hoang (2007) showed 
that the RTs amplitude may change substantially 
with the angle between the radiation direction and 
magnetic field. 
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Fig. 4. Grain alignment in surface layer. The ratio 
(^rad/^T) 2 \^ a exceeds 10 when r > 1AU, which means 
that some grains in the surface layer are aligned when 
r > 1AU. Results are for a — 1/im grains. From Cho & 
Lazarian (2007). 



Fig. 5. Illustration of grain alignment on disk surface. 
Near the star, grains are not aligned due to high density. 

4. GRAIN ALIGNMENT IN DISKS 

We use Eq. ©, instead of the DDSCAT software 
package, to obtain radiative torque (Qr) on grain 
particles in the T Tauri disks. We take a conservative 
value of Qr at A ~ a: Qr ~ 0.1 at A ~ a. Apart 
from Qy, we also need to know u ra( i and to get 
the (w ra d/wT)La ratio (see Eq. ([5])). We directly 
calculate u ra d and nji using the disk model in C01. 
We assume that T drag - T dragsas . 
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We assume that the RT alignment is perfect when 
the ratio (u r ad/wT)x~a exceeds lfJE Calculations 
(see Cho <fc Lazarian 2007) show that grains near 
the central star cannot be aligned due to high gas 
density near the star. Indeed Figure 2] shows the ra- 
tio (oJrad/ujT^r^a exceeds 10 when the distance from 
the central star, r, is large, which means that grains 
in the surface layer are aligned when r is large. We 
expect that polarized emission from the surface layer 
is originated from outer part of the disk. Similar re- 
sults hold true for disk interior. Calculations show 
that, at large r, large grains are aligned even deep 
inside the interior. On the other hand, at small r, 
only grains near the disk surface are aligned (see an 
illustration in Figure [5]). 




1 10 100 1000 

Wavelength (/jm) 



Fig. 6. Spectral energy distribution. The vertical axis 
(i.e. \F\) is in arbitrary unit. Thick solid line: total 
(i.e. interior + surface) emission from disk. Thin dotted 
line: total emission from disk surface. Thick dotted line: 
polarized emission from disk surface. Thin dashed line: 
total emission from disk interior. Thick dashed line: po- 
larized emission from disk interior. Note that, in these 
calculations of polarized emission, we ignored the direc- 
tion of polarization vectors and we only take the absolute 
value of them. Results for oblate spheroid grains with 
axis ratio of 1.3:1. From Cho & Lazarian (2007). 



3 The perfect alignment is true for grains having superpara- 
magnetic inclusions (Lazarian & Hoang 2008). For ordinary 
paramagnetic grains the degree of alignment can vary. 



5. PREDICTIONS FOR DEGREE OF 
POLARIZATION 

5.1. Estimates for Spectral Energy Distribution 

In this subsection, we calculate the degree of po- 
larization of emitted infrared radiation from a disk 
with structure and parameters described in C01. In 
this subsection, we assume that the disk in face-on. 
The degree of polarization will be zero for a face-on 
disk when magnetic field is perfectly azimuthal and 
the disk is cylindrically symmetric. In this section, 
we are concerned only with the absolute magnitude 
of the polarization. 

Figure [5] shows the results for 1.3:1 oblate 
spheroid. The degree of polarization can be as large 
as ^5% in FIR/sub- millimeter wavelengths and ~ 
10% in mid-IR regimes. The polarized emission at 
FIR is dominated by the disk interior and that at 
mid-IR is dominated by the disk surface layer. Note 
again that, in these calculations, wc ignored the di- 
rection of polarization and we only take the absolute 
value of it. Note that, since the degree of polariza- 
tion of emission from the disk surface layer is very 
sensitive to the maximum grain size in the surface 
layer, the results for A < 100/xm should be very sen- 
sitive to the maximum grain size in the surface layer. 

5.2. Radial energy distribution 

Figure [7] shows radial distribution of emitted ra- 
diation. For A = 850/xm, both radiations from the 
disk interior and the surface layer are dominated by 
the outer part of the disk. But, for A = 10/im, the 
inner part of the disk contributes significant portion 
of total emission and polarized emission. 

5.3. Effects of disk inclination 

In this subsection we calculate actual degree of 
polarization that we can observe. Chiang & Gol- 
drcich (1999) calculated spectral energy distribu- 
tion (SED) from inclined disks. We follow a similar 
method to calculate the the SED of polarized emis- 
sion. 

Figure [8] shows the effects of the disk inclination. 
We calculate the polarized emission from the disk 
interior. The viewing angle 9 (=the angle of disk 
inclination) is the angle between the disk symmetry 
axis and the line of sight. We plot the direction of 
polarization for 3 different wavelengths and 2 differ- 
ent viewing angles. The lines represent the direction 
of polarization. Since we assume that magnetic field 
is azimuthal, the direction of polarization is predom- 
inantly radial (see lower panels). 
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Fig. 7. Radial energy distribution, (a) A = 10/j,m. Inner part of the disk emits substantial amount of radiation. But 
it emits negligible amount of polarized radiation. Note that, when r < 1AU, grains in the surface layer are not aligned 
and only negligible fraction of grains are aligned in the interior (see Figure 2] for disk surface), (b) A = 50/xm. (c) 
A = 100/xm. (d) A = 850/im. The result for A = 450/im (not shown) is very similar to that for A = 850/im. From Cho 
& Lazarian (2007). 



6. PROSPECTS 

Multifrequency observations of protostellar disks 
have become a booming field recently. As grains at 
different optical depths have different temperatures, 
multifrequency measurements reveal the structure 
of the disk. They have advanced substantially our 
knowledge of the disks and allowed theoretical ex- 
pectations to be tested. 

Our study reveals that multifrequency polarime- 
try is very important for the protostellar disks. 
The synthetic observations that we provide explic- 
itly show that observations at wavelength less than 
100 (im mostly test magnetic fields of the skin lay- 
ers, while at longer wavelengths test magnetic fields 
of the bulk of the disk. Therefore polarimetry can, 
for instance, test models of accretion, e.g. layered ac- 
cretion (Gammie 1996). Combining the far- infrared 
polarimetry with polarimetric measurements at dif- 
ferent frequencies may provide additional insight 



into the magnetic properties of protostellar accretion 
disks. 

Most of the present day polarimetry will be done 
for not resolved protostellar disks. The size of the T 
Tauri disks is usually less than ~ 300 AU (see, for 
example, C01). If we take the distance to proto-stars 
to be around > lOOpc, then the angular sizes of the 
disks are usually smaller than 6". The angular reso- 
lution of SCUBA polarimeter (SCUPOL) is around 
14" (Greaves et al. 2000) and that of SHARC II po- 
larimeter (SHARP; Novak et al. 2004) at 350/im is 
around 9". Therefore it is not easy to obtain plots 
like Figures[H The angular resolution of the intended 
SOFIA polarimeter is around 5" at 53fim, 9" at 
88/ito, and 22" at 215/im. We see that the intended 
SOFIA polarimeter will be at the edge of resolving 
structure of close- by disks, while other instruments 
will not resolve typical T Tauri disk. Therefore for 
most of the near future observations our predictions 
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Fig. 8. Simulated observations. Degree of polarization is calculated for the total radiation (i.e. interior + surface) 
from the disk. The disk inclination angle 6 is the angle between disk symmetry axis and the line of sight. From Cho & 
Lazarian (2007). 



in Fig. [5] and [5] are most relevant. 

7. SUMMARY 
Making use of the recent advances in grain align- 
ment theory we calculated grain alignment by RTs in 
a magnetized T Tauri disk. Based on this, we calcu- 
lated polarized emission from the disk. Our results 
show that 

• Polarization arising from aligned grains reveals 
magnetic fields of the T Tauri disk. 

• Disk interior dominates polarized emission in 
FIR/sub- millimeter wavelengths. 

• Disk surface layer dominates polarized emission 
in mid-IR wavelengths. The degree of polariza- 
tion is very sensitive to the maximum size of 
grain in the disk surface layer. 

• Our study of the effect of the disk inclination 
predicts substantial changes of the degree of po- 
larization with the viewing angle. The coming 



mid-IR/FIR polarimeters are very promising for 
studies of magnetic fields in protostellar disks. 

• Polarization at different wavelengths reveals 
aligned grains at different optical depths, which 
allows one to tomography magnetic field struc- 
ture. 
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